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Abstract 

We give examples where the Heun function exists as solutions of wave 
equations encountered in general relativity. As a new example we find 
that while the Dirac equation written in the background of Nutku 
helicoid metric yields Mathieu functions as its solutions in four space- 
time dimensions, the trivial generalization to five dimensions results 
in the double confluent Heun function. We reduce this solution to the 
Mathieu function with some transformations. 
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1 Introduction 



Most of the theoretical physics known today is described by a rather few num- 
ber of differential equations. If we study only linear problems, the different 
special cases of the hypergeometric or the confluent hypergeometric equa- 
tion often suffice to analyze scores of different phenomena. These are two 
equations of the Fuchsian type with three regular singular points and one 
regular, one irregular singular point respectively. Both of these equations 
have simple recursion relations between two consecutive coefficients when a 
power expansion solution is attempted. This fact gives in many instances 
sufficient information on the general behavior of the solution. If the problem 
is nonlinear, one can usually fit the equation describing the process to one of 
the different forms of the Painleve equations [lj. 

It seems to be a miracle that so diverse phenomena, examples in potential 
theory or wave equations with physical applications, can be described with so 
few equations. Physicists are lucky since most of the phenomena in physics 
of the present can be described in terms of these rather simple functions. 
Perhaps, this refers to a symmetry beyond all these things, like the occur- 
rence of hypergeometric functions may signal to the presence of conformal 
symmetry. 

In the linear case, sometimes it is necessary to go to equations with more 
singular points. The Heun equation [2] [3] |1] , its confluent cases, or its special 
cases, Mathieu, Lame, Coulomb spheroidal equations etc., all have additional 
singular points, either four regular, or two regular, one irregular, or two ir- 
regular. The price you pay is the fact that in general there exists no recursion 
relation between two consecutive coefficients when a power series expansion 
is used for the solution. A three or four way recursion relation is often awk- 
ward, and it is not easy to deduce valuable analytical information from such 
an expansion. New versions of the computer package Maple, for example 
Maple 10, gives graphical representations of Heun functions. Although this 
is a great help, it is sometimes awkward to use these functions as potentials 
in either wave or Schrodinger equations. This may be a reason why much 
less is known about these equations compared to hypergeometric functions 
and all the other functions derived from them. 

One encounters Mathieu functions when one uses elliptic coordinates, in- 
stead of circular ones, even in two dimensions [5|. Phenomena described by 
Heun equations are not uncommon when one studies problems in atomic 
physics with certain potentials [EJwhich combine different inverse powers 
starting from the first up to fourth power or combining the quadratic po- 
tential with inverse powers of two, four and six, etc. . They also arise 
when one studies symmetric double Morse potentials. Slavyanov and Lay [7j 
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describe different physical applications of these equations. Atomic physics 
problems like separated double wells, Stark effect [H],[Zj, hydrogen- molecule 
ion [9], [7j use these forms of these equations. Many different problems in 
solid state physics like dislocation movement in crystalline materials, quan- 
tum diffusion of kinks along dislocations are also solved in terms of these 
functions [7J. The famous Hill equation [TU] , [T] for lunar perigee can be cited 
for an early application in celestial mechanics. 

In general relativity, while solving wave equations, we also encounter dif- 
ferent forms of the Heun equations . Teukolsky [H],[Tj studied the pertur- 
bations of the Kerr metric and found out that they were described by two 
coupled singly confluent Heun equations. Quasi-normal modes of rotational 
gravitational singularities were also studied by solving this system of equa- 
tions [12], [Ij. 

In recent applications they become indispensable when one studies phe- 
nomena in higher dimensions, for example the article by G. Siopsis [13] . or 
phenomena using different geometries. An example of the latter case is seen 
in the example of wave equations written in the background of these met- 
rics. For instance, in four dimensions, we may write wave equations in the 
background of 4D Euclidean gravity solutions. For the metric written in 
the Eguchi-Hanson instanton [H] background, the hypergeometric function 
is sufficient to describe the spinor field solutions [TH1HE]- One, however, has 
to use Mathieu functions to describe even the scalar field in the background 
of the Nutku helicoid instanton (TTJ HH] when the separation of variables 
method is used for the solution. Schmid et al[19j have written a short note 
describing the occurrence of these equations in general relativity. Their ex- 
amples are the Dirac equation in the Kerr-Newman metric and static per- 
turbations of the non extremal Reisner- Nordstrom solution. They encounter 
the Generalised Heun Equation [3J[4J[20J while looking for the solutions in 
these metrics. Here we see that as the metric becomes more complicated, one 
has to solve equations with larger number of singular points, with no simple 
recursion relations if one attempts a series type solution. As a particular 
case of confluent Heun equation, Fiziev studied the exact solutions of the 
Regge- Wheeler equation |21j[22j. One also sees that if one studies similar 
phenomena in higher dimensions, unless the metric is a product of simple 
ones, one has higher chances of encountering Heun type equations as in the 
references given [23] [24] . 

Here we want to give further examples to this general behaviour. Our 
first example will be the case already studied by Sucu and Unal [T5] . In their 
paper they study the spinor field in the background of the Nutku helicoid 
instanton. They obtain an exact solution, which, however, can be expanded 
in terms of Mathieu functions [25]. At this point note that their solution 
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in the background of the Eguchi-Hanson instanton is expanded in terms of 
hypergeometric functions. Taking the next 4D Euclidean gravity solution in 
row, the Nutku helicoid instanton with a two centers [26] [27], results in a 
function with higher singularity structure. The second example is the similar 
equation in five dimensions. Here the metric is extended to five dimensions 
with a simple addition of the time coordinate. In this case we obtain the 
double confluent Heun function as the solution, which can be reduced to the 
Mathieu function with a variable transformation. 

At this point the motive for our paper becomes clear. Both in four and 
five dimensions, in the background metric of the Nutku helicoid instanton, we 
can write the solutions of the Dirac equations in terms of Mathieu functions. 
There is a catch here, though. Although we can express the solution in a 
closed form in the four dimensional CclSG, clS done by Sucu and Unal [To] , 
this is not possible using the solutions in five dimensions, since the constants 
used in the two equations are not the same. 

In four dimensions we can also calculate the the Greens function for this 
differential equation in closed form following the steps in reference [18]. In 
five dimensions we could not succeed in summing the the infinite series of the 
product of Mathieu functions to express both the solution to the differential 
equation and the Greens function of the same differential equation in a closed 
form. This is due to the existence of two different constants in the two 
Mathieu functions used in the expansion. What we wanted to show is that 
going to one higher dimension, we got solutions which were more complicated. 

We describe our examples in the consecutive sections, first the case in 
four, then in five dimensions. We then give the solutions for the scalar 
operator. We end with some additional remarks. In our work we use only 
the massless field, since taking the massive field is technically like going one 
higher dimension, which complicates the problem. We instead go to one 
higher dimension in an explicit fashion in the consecutive section. 

2 Equations in Four Dimensions 

The Nutku helicoid metric is given as 



ds 2 



1 



[ rfr 2 + ( r 2 + a 2 )M 2 + 1 + ^ S in 2 9 dy 2 
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where < r < oo, < 8 < 2ir, y and z are along the Killing directions and 
will be taken to be periodic coordinates on a 2-torus |18j . This is an example 
of a multi-center metric. This metric reduces to the flat metric if we take 
a = 0. 

ds 2 = dr 2 + r 2 d9 2 + dy 2 + dz 2 . (2) 
If we make the following transformation 

r = a sinh x, (3) 

the metric is written as 

a 2 

ds 2 = — sinh 2x(dx 2 + d6 2 ) 
2 

H [(sinh 2 x + sin 2 9)dy 2 (4) 

sinh 2x 

— sin 29dydz + (sinh 2 x + cos 2 0)dz 2 ] . 



We use the NP formalism [2EJ EHj in four Euclidean dimensions 
To write the Dirac equation in this formalism, we need the choose the base 
vectors and calculate the spin coefficients, the differential operators and the 
7 matrices in curved space. 
We take base vectors 

e£ = {P,P,m^,m^}, (5) 



to give 
The tetrad 
satisfies 



ds 2 = l®l + l®l + m®m + m®m. (6) 
e a = e a v dx\ (7) 



where r] a b is the flat metric. 
We choose 

r = ^=L=(l,z,0,0), (9) 

ra* 1 = (0,0,cosh(x - i6), isinh(x - i6)), (10) 

V sinh 2x 

giving the two non zero spin coefficients 

cosh (2a;) 



a sinh 3/2 2x ' 



'ID 
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a sinh 3/2 2x ' 



The rest of the spin coefficients, 

K = v = ~f = a = j3 = 'K = T = ii = \ = p = 0. 
These expressions give the differential operators 



D = rrfd,, = 



1 



D = fn^d. 



a/ sinh 2x 
1 

V sinh 2x 



5 = m 



[cosh(x — id)d y + i sinh(x — i9)d z \ , 

[cosh(z + i8)d y — i sinh(z + i9)d z \ , 
[d x + id e ], 



1 



5 = Fd, 



aV sinh 2x 
1 



^ aV sinh2x 
The massive Dirac equation reads 



[d x - ide] 



where 



The 7 matrices can be written in terms of base vectors as 
7" = V2 



( V 

-rh^ P 







V 



rrv 



n in 



/ 



These matrices satisfy 

{-f,-f} = 2cr- 

The spin connection is written as 



In expanded form, these equations read 



V2 



aVsinh2x 



{(d x + id e )* 3 



r m ■ no • //i ■ no it MaV sinh 2x T , 

+a[cos(# + lajdj, + sm(6> + tx)d z ]^ A -= ^} = 0, 

V2 
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V2 {iOr • 



a\J sinh 2x 



r ,„ . ,„ . . _ . _ Ma V sinh 2x T , „ 

-a[cos(# - za;)d y + sin(0 - n;)<9 2 ]^ 3 — = ^ 2 } = 0, (24) 

V2 



a/2 

{(a x - ia e +coth2z)^ 1 



aV sinh 2x 



r . s« . //, . x^ lT May/ sinh 2x T , „ 

-a[cos(# + zx)^ + sin(0 + ix)d z }^ 2 7= ^3} = 0, (25) 

V2 



/2 

(9 X + i<9 e + coth2x)# 2 



a\/ sinh 2a; 

+a[cos(# - + sin(# - 1 - Mav/s *? h 2x ^ 4 = . (26) 

v2 

To simplify calculations, we will study the massless case. Then we see that 
only {^1, ^2} and {^3, ^4} are coupled to each other. If we take 

^ = e *(M+*.*)^ i ( x> 0) (27) 
and make the transformations 

k y = k cos (f),k z = k sine/), (28) 



we get 



and 



Now we have to solve 



sinh [a; — i(6 — 0)1 . _ 

*i = r^r < 29 ) 

V sinh 2x 

sinh fx + i(0 — 6)] T . . 

* 2 = L . V Vn ^> 2 . (30) 

V sinh 2a; 



-2 f a 2 A; 2 1 

£,,2*1,2 = 2afc ^^ + ^ + — {cos[2(^ + 0)] - cosh 2a;} j * 1>2 = 0, 

(31) 

whose solutions can be expressed in terms of Mathieu functions. 



7 



*i = e 



lk(z sin 4>+y cos </>) Smh[z j{ 



V sinh 2x 



a 2 k 2 



x 



Sefa,- 
5e(Ci 



■x 



a 2 A: 2 



4 

a 2 A; 2 



+ 5o(Ci, — 



21.2 



2 7 



h 



(32) 



*5 



= jfc(2 sin 4>+y cos < 



, sinh[x + i{9 — <fi)] 



V sinh 2 



a 2 A; 2 

Se(C 2 , — ,-ix) + So(C 2 , 

a 2 A; 2 



x 

a 2 A; 2 



4 



a 2 A; 2 



5e(C 2 , —,0 + <t>) + Sofa, — , 9 + 0) 



}• 



(33) 



When similar transformations are done for the other components we get: 



cosh[x — i{9 
ak 



{ a 2 k 2 1 

<^ d xx + d ee + — {cos[2(# + 0)] - cosh 2x} \ ^ 3 = 

(34) 



cosh[x + i{9 — 0)1 f. n a 2 k 2 r , , n , ,1 
^7 — | 4* + + — {cos[2(# + 0)] - cosh2:r} \ ^ 4 = 

(35) 



The solutions again can be expressed in terms of Mathieu functions: 



■qf^ _ e ifc(2sin0+?/cos0)| 



a 2 A; 2 a 2 k 2 
Sefa, -is) + So(C 3 , — : 



x 



^(C3,-^^ + 0) + 5o(C3 



a 2 A; 2 



;0 + <f>) 



(36) 



3 ifc(zsin0+ycos< 



a 2 A; 2 a 2 A; 2 
Sefa, — , —ix) + Sofa, — , —ix) 



x 



n 2 k 2 n 2 k 2 

Sefa, 9 + 0) + Sofa, — , 9 + 0) 



4 
}• 



(37) 
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Here note that 



a 2 k 2 



Se(C 



4- 



-ix) + So((. 



2 1,2 



a 2 k 



4 ; 



-ix) 



can be expressed in terms of modified Mathieu functions with real arguments. 
Here Q are separation constants. 

At this point also note that we can get solutions for ^4 in the plane wave 
form, which are given as exp(ka(sin(9 — <fi + ix) + sin(9 — <fi — ix))), similar 
to the ones given by Sucu and Unal [IS]. Here, since we want to point to the 
occurrence of Mathieu functions in mathematical physics, we use the product 
form. This form is also more useful when boundary conditions are imposed 
on the solution. 



3 Equations in Five Dimensions: 



The addition of the time component to the previous metric gives: 

ds 2 = —dt 2 + ds 2 , 
resulting in the massless Dirac equation as: 



Here 
and 



(7^ + 7^-7^-7*1^ = 0. 
T t = 



7 



/ i \ 

i 

-i 

\ -i J 



(38) 

(39) 
(40) 

(41) 



giving the set of equations, 



y/2 



iy/ sinh 2 



-.{{d x + id e )% 



X 



cl\I sinli 2x 

+a[cos(# + ix)d y + sin(6» + ix)d z }^± + i -= d t ^i} = 0, (42) 

v2 
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y/2 



ay/ sinh 2x 



{(d x -id e )V 4 



-a[cos(6 - ix)d y + sm{9 - ix)d z ]^ 3 + i ay/sin ^ 2x d t ^ 2 ] = 0, (43) 

v2 



y/2 



a\J sinh 2x 



{(d x - id e + coth2x)^i 



Qj\i smli 2x 

-a[cos(9 + ix)d y + sin(6> + ix)d z }^> 2 - i 7= d t ^ 3 } = 0, (44) 



V2 



a\J sinh 2x 



{(d x + id e + coth2x)^ 2 



+a[cos(9 - ix)d y + sin(6> - ix)d z ]^ 1 - z -= <9 t * 4 } = 0. 

a/2 



(45) 



If we solve for ^ and ^2 and replace them in the latter equations, we get 
two equations which has only ty 3 and ^4 in them. If we take 



the resulting equations read: 



(46) 



a 2 k 2 

dxx + d 8e + — ^— {cos[2(6> + </>)] - cosh 2x} + 2a 2 k 2 sinh 2x \ ^ 3i4 = 0. 

(47) 

If we assume that the result is expressed in the product form ^3 = Ti(x)T 2 (9), 
the angular part is again expressible in terms of Mathieu functions. 



T 2 {9) = Se 
+So 



a 2 k 2 . /l + cos(fl + 0). 

7], — , arccos( 1 



V, 



4 
a 2 k 2 



,l + cos(# + 0). 
, arccos(W ) 



(48) 



Here i] is the seperation constant and the periodicity on the solution makes 
it equal to the square of an integer. 
The equation for T\ reads: 
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da k 9 7 9 



cosh 2x + 2a fc t sinh 2x — ,7 > Ti = 



(49) 



The solution of this equation is expressed in terms of double confluent Heun 
functions [I]. 



71 Or) 



HeunD 
+HeunD 



a k ,o & k 

0, — — h 77,4a r/c; 



2 

a 2 k 2 

0, 77 + ^, 4a 2 A; 2 ,— 
—dx 



a 2 _k 2 

y 



— 77, tanh a; 
77, tanhx 



(50) 



1 2 • 



HeunD [0,rj + ^,4a 2 A; 2 , ^| 77, tanh x] 

We only take the first function and discard the second solution. We see 
that as x goes to infinity, the function given above diverges. The function 
is finite at x = though. In order to get well defined functions, we study 
the region where x < F, where F is a finite value. We will give a way to 
determine F below. 

We can use either Dirichlet or Neumann boundary conditions for our problem 
in four dimensions. There is an obstruction in odd Euclidean dimensions that 
makes us use the Atiyah-Patodi-Singer [33] spectral boundary conditions. 
These boundary conditions can also be used in even Euclidean dimensions if 
we want to respect the charge conjugation and the 7 s symmetry [34] . 
Just to show the differences with the four dimensional solution, we attempt 
to write this expression in terms of Mathieu functions. This can be done 
after few transformations. We define 



A = 2a z k 



2?„2 

t > 



(51) 



and use the transformation 



B = - V , 



C 



a 2 k 2 



2 ' 



z = e 



-2x 



(52) 
(53) 

(54) 



Then the differential operator is expressed as 
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= Az 2 d zz + Azd z + A'z + B + C 



Here 



A 



, C- A 



C 



, C + A 



(55) 

(56) 
(57) 

(58) 
(59) 

2-4- ™ 
The solution of this equation is also expressible in terms of Mathieu functions 
given as: 



If we take 



and 




u — z, 



w = -(u+-) 



and set E = y/A'C we get, 



E B 

= {w 2 - l)d ww + wd w + — w + — . 



R(z) = Se(—B, E, arccos s/w + 1) + So(—B, E, arccos yw + 1) (61) 

At this point we see a natural limitation in the values that can be taken by 
our radial variable since the argument of the function arccos can not exceed 
unity. The fact that + 1 can not exceed unity limits the values our initial 
variable x can take, thus determining F which imposed on our solution in 
equation [50]. 

Here we also see a difference from the four dimensional case. Although both 
the radial and the angular part can be written in terms of Mathieu functions, 
the constants are different, modified by the presence of the new —2a 2 k 2 2 term, 
which makes the summation of these functions to form the propagator quite 
difficult. 

In four dimensions we can use the summation formula [351 EH] for the product 
of four Mathieu functions, two of them for the angular and the other two for 
the radial part, summing them to give us a Bessel type expression. This result 
makes the calculation of the propagator, similar to the case given in reference 
[18] possible. In that case, the similar analysis also gives the solution to the 
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differential equation in a closed form as given in reference [15]. Here since 
the radial and the angular parts have different constants, this summation 
formula is not applicable write the Greens function in a closed form. We also 
see that using the generating function formula for these functions to write the 
solution to the differential equation in terms of plane waves, as described in 
the paper by L.Chaos-Cador et al., is not applicable in the five dimensional 
case due to the same reason. 



4 Laplacian 

In this section we give the Laplacian operator written in this background. It 
is used for the calculation of the field equation for the scalar particle, similar 
to the case studied in reference [18]. 

H : = -L&yz^fy (62) 

H : = -d tt + d xx + d ee + a 2 sinh 2 2x(d yy + d zz ) 

+a 2 (cos 0d y + sin 6d z ) 2 — a 2 sinh x cosh xd u (63) 

We see that there are three Killing vectors and one quadratic Killing tensor 
with eigenvalues given below. From the Killing tensor we can construct a 
second order operator tensor [18] 

K = ~d ee ~ a 2 (cos 9d y + sin 9d z f (64) 

K§ = A$. (65) 
We use A as the separation constant. The other eigenvalues are 

<9 t $ = kt$, (66) 

d y $ = k y $, (67) 
d z § = k z §. (68) 

We have 

$ = e i{ktt+k y y+k * z) R(x)S(6) (69) 
where S obeys the equation 

d 2 S 



dO 



^ + (\-a 2 k 2 cos 2 8)S = 0. (70) 
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Here 

~e = e-(f) (71) 

and the solution reads 

S(6>) = Se(^ + A,^f ,0-0) 

+ 5o(-^- + A,^L*-0. (72) 

A is the separation constant which goes to the square of an integer due to 
periodicity of the function. 
The radial part obeys 

d 2 R A 

—-— + a 2 (sinh x cosh xkf — k 2 sinh 2 x -)R = 0. (73) 

dx 2 a 2 

The solution to this equation, the radial solution, can be written in terms of 
Double confluent Heun functions, 



R(x) = HeunD [0, a 2 k 2 — A, —a 2 k 2 , A, tanhx] 

+HeunD [0,a 2 k 2 - A, -a 2 k 2 , A,tanh x] (74) 

/—dx 
9 2- 
HeunD [0, a 2 A; 2 — A, —a 2 k 2 , A, tanh x] 

and this can be reduced to the modified Mathieu function after performing 
similar transformation as in the spinor case treated above. Taking A = 
B = — A and C = ~ a ^ fc2 , and use z = e~ 2x transformation we get the 
same result as in equation (61). 



5 Conclusion 

Here we related solutions of the Dirac equation in the background of the 
Nutku helicoid solution in five dimensions to the Double confluent Heun 
function. This solution can be also expressed in terms of the Mathieu func- 
tion, which is more familiar to the physics community, at the expense of 
using the z = e~ 2x transformation, which maps infinity to zero followed by a 
rescaling and a further transformation where, aside from the scaling, we are 
taking the hyperbolic cosine of the original radial variable. This transforma- 
tion also brings a natural limit on the radial variable. Essentially we do not 
gain much, since the Mathieu function also has a three way recursion rela- 
tion, and a not very handy generating function, compared to the generating 
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functions belonging to the hypergeometric family. Such a transformation to 
Mathieu functions may not be possible when more complicated backgrounds 
are taken. Although in many cases the Heun function may reduce to more 
simple forms [36] [37], there are metrics in general relativity where this may 
not be possible. The Nutku helicoid metric is such an example. 
We know that using Maple 10 programme we can see the graphical represen- 
tation of these functions. Different asymptotic solutions are also studied in 
references [3,4,5] and in individual papers. Still it is our feeling that most of 
the mathematical physics community is not at ease with these functions as 
they are with the better known hypergeometric or confluent hypergeometric 
functions. The literature on these functions is also very limited. We think we 
have covered most of the monographs where these functions are thoroughly 
studied in our bibliography. If we search the SPIRES web site, which is 
commonly used by mathematical physicists, the number of entities is rather 
small. This is the reason we think one should be exposed to its applications 
more often. 

We also wanted to show in our work that when one uses more complicated 
forms of similar structures, one gets more complicated solutions. What we 
mean by this phrase that if one writes the Dirac equation in the background 
of the simplest Euclidean solution in four dimensions, the Eguchi-Hanson 
solution, one gets the hypergeometric function as solutions of the wave equa- 
tion [TSj. If we use the next solution in the order of complexity, the Nutku 
helicoid solution, one gets a solution of the Heun type. We think the solu- 
tion obtained in reference [15] for this case in terms of exponential functions 
is somewhat misleading, since there they get the generating function of the 
Mathieu functions. This may not be recognized by the readers unless one is 
an expert in this field. 

We also studied the same equation in one higher dimension. Often increasing 
the number of dimensions of the manifold in which the wave equation is 
written, results in higher functions as solutions. Here we call a function of a 
higher type if it has more singularities. In this respect Heun function belongs 
to a higher form than the hypergeometric function. In the Introduction we 
gave examples of the use of Heun functions encountered in different physics 
problems. 

In the future we will try to find further examples of such functions encoun- 
tered in solutions to wave equations in general relativity. 
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